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SUMMARY 


A study of single-degree-of -freedom pitching oscillations of a wing 
has been presented. This study includes the effects of Mach number and 
struct\u*al damping and is primarily an extension of a recent paper by 
Smilg* in which incompressible flow was considered. The actual existence 
of single-degree-of -freedom flutter was demonstrated by some low-speed 
tests of a wing, pivoted a short distance ahead of the leading edge with 
a geometric aspect ratio of 5»8?. In general, good agreement was found 
between experimental and calculated results for high values of an inertia 
parameter corresponding to high altitudes, but differences exist for low 
values of the inertia parameter. The effect of aspect ratio has not 
been considered in the calculations and could have an appreciable 
influence on the oscillation. 


INTRODUCTION 


The possibility of the existence of single-degree-of -freedom 
oscillatory instability or flutter in incompressible flow, both potential 
and separated, has been known for some time. As early as 1929 Glauert 
(reference 1) noted the possible loss of damping of a pitching w in g in 
incompressible flow which might lead to an oscillatory instability that 
may be referred to as single-degree-of -freedom flutter. In 1937 Possio 
made similar observations for supersonic flow (reference 2) and in 19 kS 
this study was elaborated on by Garrick and Rubinow (reference 3), who 
observed that, under certain conditions, a single-degree- of -freedom 
oscillation is possible in incompressible flow. Subsequently, Smilg 
(reference U) made calculations showing the ranges of axis-of -rotation 
location and an inertia parameter which could lead to an oscillatory 
instability in pitch or yaw for the incompressible case. 

Until recently whatever interest was shown in single-degree-of- 
freedom flutter was largely academic because the ranges of parameters 
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involved did not appear practical. However, with current airplanes and 
missiles designed for high speeds and high altitudes, the subject becomes 
a more practical one, for under these conditions undamped oscillations 
of even very small amplitude may become important. In addition, calcu- 
lations of single-degree-of -freedom flutter may represent a useful, 
easily obtained limit for cases of coupled flutter involving other 
degrees of freedom. 

This paper considers specifically the type of single-degree-of- 
freedom flutter associated with the pitching of an airfoil about various 
locations of the axis of rotation. It extends the work of reference U 
to include the effects of Mach number up to M = 0.7 and discusses the 
effect of structural damping for one location of the axis of rotation. 

The results of an experimental investigation which confirms the existence 
of single-degree-of -freedom flutter are compared with the theoretical 
values. The calculations were based on two-dimensional aerodynamic- 
force coefficients and involved a single degree of freedom. The effect 
of aspect ratio and the coexistence of other degrees of freedom would 
modify the resiilts to a large extent. 


SYMBOLS 


a nondimensional distance from midchord to axis of 

rotation, based on half-chord, positive rearward 

b half-chord 

c spring constant 

Ca coefficient of torsional rigidity per unit length 

d damping coefficient 

F and G functions of k for oscillating plane flow 

ga structural damping coefficient 

la moment of inertia about axis of rotation per unit 

length 

laa out-of -phase (imaginary) component of moment on 

airfoil about axis of rotation per unit length 


k 


reduced frequency (bco/v) 
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m 

mass 

M 

Mach number 

Ma 

aerodynamic moment per unit length 


real part of aerodynamic moment, for incompressible and 
compressible flow, respectively 


aerodynamic flutter coefficients (see reference 5) 

^aa 

in-phase (real) component of moment on airfoil about 
axis of rotation per unit length 

V ■ 

flutter velocity 

X 

displacement 

a 

angrilar displacement about axis of rotation, positive 
in stalling direction 

P 

fluid density 

00 

circular frequency at flutter 

■ 

natural circular frequency 


Al'IALYSIS 

Introductory Considerations 


Before the specific example of single-degree-of -freedom pitching 
flutter is discussed, it may be advantageous first to review the concept 
of a single-degree-of -freedom vibrating system and then to show the 
relation of this example to an aerodynamic system. The linear differen- 
tial equation for a free system consisting of a mass m,^a spring having 
a spring constant c, and a viscous damper having a coefficient d is 

mx + dx + cx = o n'i 


The motion represented by this equation is damped if d is a 
positive quantity, as is ordinarily the case. If d should be negative. 
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the motion is undamped, a condition of dynamic instability exists, and, 
if d is zero, harmonic oscillations corresponding to a borderline 
condition between damped and undamped motion may exist. 

For a system such as an aircraft wing oscillating in a steady air 
stream, the same type of equation would apply as for the mass-spring- 
damper system previously mentioned. However, the coefficients m, d, 
and c of equation (1) will now have added components associated with 
the aerodynamics. The equation for a wing oscillating in pitch in a 
steady two-dimensional air stream is 


laH + (1 + ig(j)Caa = Ma(a,d,a, . . . ) 


( 2 ) 


where represents the complex aerodynamic moment, which is a function 

in part of amplitude a, velocity d, acceleration a, reduced frequency 

k = location of axis of rotation a, Mach number M, and sweep angle. 

Equation (2) is complex and may be separated into two components, one 
associated with the damping of the system (sometimes called the imaginary 
part) and the other associated with the flutter frequency and velocity 
(sometimes called the real part). 


Equation for Pitching Oscillations, M = 0 

From reference 6, the values of the damping equation and the 
frequency equation for two-dimensional incompressible flow are: Damping 

equation: 



Equation (3) is equivalent to the vanishing of the damping coefficient 
d of equation (1) and thus represents a borderline condition between 
damped and undamped oscillations. The flutter frequency and velocity 
may then be determined from equation (It) . 
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Equation (3) cannot be solved explicitly for k since the functions 
F and G are transcendental functions of k. There are several methods 
of solving this equation; a convenient one, given in reference 7, is to 
assume values of l/k and solve for the structural damping coefficient. 
The type of structural damping force commonly used in flutter calculation 
is in phase with the velocity but. proportional to the amplitude. If the 
damping coefficient is plotted against l/k, the value of l/k for aity 
given damping coefficient may be determined. When the value of l/k 
that satisfies the imaginary or damping part of the moment equation has 
been determined, the frequency of oscillation and the velocity may be 
determined from the real part of the moment equation, equation (U). 
Equation (i;) may be put in different form as follows: 



where 

Mr 


If the torsional restraint coq is zero, equation (5) reduces to 



I 

np 



r 


( 6 ) 


so that, if the value of the inertia parameter Ig/nph^ exceeds the 
value of Mr for the given axis-of -rotation location, the oscillation 
can exist at all airspeeds above zero speed. The frequency is then a 
direct function of the velocity as defined by the following equation: 


V = 


to 

k 


(7) 


where l/k is the value of the flutter-speed parameter associated with 
the borderline condition between damped and undamped oscillation for 
the given axis of rotation. 


Equation for Pitching Oscillation Including Mach Number 

In order to consider the effect of Mach number, the results of 
reference 5 may be used. The method of computation is the same as 
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described in the preceding section; however, the aerodynamic moment Mq 
has been redefined to include the effect of Mach number. 


The damping (imaginary) con^jonent (see references 5 and 6 ) is 


i-aa 




and the frequency (real) equation is 


.^2 = 0 


( 8 ) 


CO _ 

“a/ 


1 - Mr' 


Tipb^ 


where 


Mr' 


k^L 


”3 - P 


11 


2\2 


+ a 


Ml + Z3 - -Zij 


1 

^1 


IJ k 


The aerodynamic coefficients Mi, M2, M3, Mj^, Zi, Z2, Z3, and Z|^ 

are functions only of reduced frequency k and Mach number (see 
reference 5). 


ANALYTICAL RESULTS 


The purpose of this section is to show the results of some calcu- 
lations made to determine the effect of some of the independent variables 
on the flutter speed and flutter frequency. 

In figure 1 , the flutter-speed parameter v/bco^j is plotted against 
the inertia parameter lo/npbU for three Mach numbers, M = 0 , M = 0 . 5 , 
and M = 0 . 7 . The region to the right and above a given curve is the 
unstable region, while the region to the left and below is the stable 
region. Increasing altitude is equivalent to increasing values of the 
inertia parameter. (Note the large change in scale between figs. 1 (a) 
to 1(f).) 

As an illustration of the meaning of the curves of figure 1 , the 
M = 0 case of figure 1 (a) (a = - 1 . 0 ) is discussed. If the value of the 
inertia parameter is below 571 (the asymptote), the configuration will 
be stable. As the altitude is increased, the inertia parameter will 
increase and, if it is equal to 571, the velocity at which an unstable 
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oscillation could occur woiold be infinite. A slight increase in the 
inertia parameter would now have a very great effect in reducing the 
critical velocity. For very large values of the inertia parameter, the 
curve is asymptotic to a value of v/bojQ which is equal to *the reduced 
velocity v/bco (that is, l/k), which for this case is 2U.72. 

The effect of Mach number is nov; examined. First, and most 
important, a large reduction in the stable region is to be noted. For 
example, in figure 1(a), the upper limit of the stable region for M = 0 
for the inertia parameter la/npbh is ^71 and this limit is reduced to 
137 at M = 0.7. 

Another effect is that, for a given speed v, the frequency of 
oscillation would increase for an increase in Mach number. For instance, 
the frequency of oscillation would be increased by a 
factor of 2 with an increase in Mach number from 0 to 0.7. 

In figure 2, the frequency ratio (co/coq)^ is plotted against the 
inertia parameter lo/npbU. This curve has the same vertical asymptote 
for the inertia parameter as for the corresponding reduced-velocity curve 
(Tig. 1) and the unstable region is again to the right and above the 
curve. The inertia parameter increases as the altitude increases. At 
low values of the inertia parameter the configuration is stable. The 
frequency of oscillation is infinite at the asymptotic value of the 
inertia parameter and decreases rapidly as the inertia parameter is 
increased further. For very large values of the inertia parameter, the 
curve is asymptotic to the natural frequency coq of the system. In 
figure 3, the minimum asymptotic value of the inertia parameter la/npb^ 
at which the oscillation could begin is plotted against Mach number for 
various positions of the axis of rotation. An important effect to be 
noted is that, as the distance of the axis of rotation from the lifting 
surface is increased, the effect of Mach number becomes increasingly 

^r63,ijiGr • 

$ 

In figure i|, the value of reduced velocity l/k is plotted against 
location of axis of rotation a for three Mach numbers. The area inside 
the curve for a given Mach number is the unstable region. The lower 
branch of each curve is asymptotic to a = -0.5 (quarter chord), but the 
upper branch has a maximum depending on the Mach number. For M = 0 
corresponding to the results of reference U, and for M = 0.5 the mLcimum 
value of a appears to be approximately -5.5; for M = 0.7 the maximum 
value of a is approximately -7. 

It should be noted that, for an airplane or missile having a 
conqjaratively short tail length (corresponding to the values of a in 
this paper), an oscillation involving a yawing motion of the vertical tail 
or a pitching motion of the horizontal tail may be an instability of the 
type considered in this paper. As a point of interest, values of the 
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inertia parameter for usual aircraft configurations when the vertical 
tail is considered as the lifting surface vary from 2,000 to 20,000 at 
sea level and would be increased by a factor of 10 for 60,000 feet. 

Since the inertia of an aircraft is usually larger about the vertical 
axis than about the horizontal axis, it appears that this type of analysis 
might be more applicable to the yawing motion. It should be noted that 
the calculations are based on two-dimensional aerodynamic coefficients 
and the effect of aspect ratio, especially if a tail surface is con- 
sidered, can be appreciable. 

In figure the effect of structural damping is shown for an axis- 
of -rotation location a = -1.2U and M = 0. The flutter- speed parameter 
v/bcoct is plotted against the inertia parameter for several values of 
structural damping coefficient g^. It is apparent that a small amount 
of structural damping has a very great effect on the flutter speed, 
especially at the low-density or high-altitude portion of the figure. 

la 

For instance, at a value of p = 18,000, a value of ga = 0.01 raises 

the flutter velocity by a factor of 3 above the zero-damping curve, and 
a value of ga = 0.02 raises the flutter velocity by a factor of 5. 
However, structura.1 damping did not influence the minimum (asymptotic) 
•value of the inertia parameter at which the oscillation could begin. 


APPARATUS AND TEST PROCEDURE 


The tests were conducted in the Langley l|.5*foot flutter research 
tunnel at low speeds (0.06<H<0.3). This tunnel can be operated at 
any pressure from atmospheric to 1/2 inch of mercury to provide a large 
range of the inertia parameter lo/npb^. 

A diagram of the model and test «section is shown in figure 6 with 
all the pertinent dimensions and parameters. The geometric aspect ratio 
was 5 . 87 . Since the wing tips were mounted close to the tunnel walls, an 
effective aspect ratio somewhat larger than the geometric aspect ratio 
was probably obtained. The wing was pivoted on ball bearings, and coil 
springs were fastened to the arms to provide structural restraint. 

A small lever was inserted through the tunnel wall and sealed with 
rubber tubing so that the wing could be disturbed while the test was 
being conducted. It was found while conducting the tests that the oscil- 
lation co\ild be started at a slightly lower fluid velocity if the wing 
were disturbed by means of the lever than if the wing were not disturbed. 
After the completion of the tests, the damping characteristics of the 
bearings were investigated, and it was found that the damping g^x was 
considerably greater for low-amplitude oscillations than for the 
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high-amplitude oscillations. This variation of damping with amplitude' 
can account for the fact previously mentioned that the model would start 
oscillating at slightly lower airspeed if the model were disturbed with 
the lever than if it were left to, the inherent air turbulence of the 
tunnel. 


EXPERIMENTAL RESULTS 


The experimental results are plotted in figure 5 where the ordinate 
is the flutter-speed coefficient v/bco^, and the abscissa is the inertia 
parameter Theoretical curves for four different values of 

damping are given, and the experimental curve is shown. 

The values of the experimental curve at the high-altitude (low- 
density) range are in close agreement with the theoretical curve for 
ga = 0,008, From an examination of the records, it appears that a 
damping coefficient of 0,01^ > ga > 0.008 was obtained} a more exact 
determination was not possible because of the dependence of the damping 
on the amplitude of oscillation. 

The important facts to be noted are, first, that a single- degree— 
of -freedom oscillation was obtained and, second, that the trend in the 
lower-density region was of the same- order of magnitude as that of the 
theoretical cujrves with damping. The reason for the discrepancy at the 
higher-density part of the plot is not known} a similar phenomenon has 
been found in other cases for the more conventional type of flutter 
involving more than one degree of freedom. 

From observations of the tests, it appears that the single-degree- 
of-freedom oscillation discussed in this paper is a mild type of flutter, 
as contrasted to the more destructive type of flutter usually associated" 
with coupled flutter. This type of instability might become of importance 
in airplane stability considerations, and the possible application to 
phenomena such as snaking should not be overlooked. It must be realized 
that three-dimensional effects may exercise some modification of these 
results. 


CONCLUSIONS 


A study of single-degree-of -freedom pitching oscillations of a wing 
has been presented. This study includes the effects of Mach number and 
structural doping and is. primarily an extension of a recent paper by 
Smilg, in which incompressible flow was considered. The actual existence 
of single-degree-of -freedom flutter was demonstrated by some low-speed 
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tests of a wing, pivoted a short distance ahead of the leading edge with 
a geometric aspect ratio of 5-87- 

The following conclusions may be drawn: 

1. The existence of single-degree-of -freedom pitching oscillations 
has been experimentally demonstrated. 

2. The experimental data are in close agreement with the theoretical 
values for high values of the inertia parameter. At low values of the 
inertia parameter, the experimental data are in poor agreement with the 
theory. 

3. Structural damping g<j has an appreciable effect on this 
instability and increases the flutter speed. 

U. The analytical results show that an increase of Mach number 
reduces the range of values of an inertia parameter for which a configu- 
ration would be stable. The results are based on two-dimensional coef- 
ficients and it is possible that aspect ratio could have a great effect. 

5 . The flutter seems to be of a mild variety, in that it would not 
necessarily cause structural failure, but the possible application to 
phenomena such as snaking for aircraft having a short tail length should 
not be overlooked. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., April 10, 199l 
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Figure 6.- Diagram of model and model installation. Inertia of system 
about axis of rotation = 0.09W foot-pound-seconds square per 
foot; natural frequency of system = 22.99. 
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